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Abstract 

The event structure associated with the production of massive photon pairs 
is studied to 0{a^as)- The cross sections for two photon and two photon 
plus one jet events are given. The DaUtz plot structure of events with two 
photons and a jet is also presented. These distributions will provide checks on 
our understanding of the production mechanisms for massive photon pairs, 
a signal thought to be important for searches for the Higgs boson in the 
intermediate mass region. 
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I. INTRODUCTION 



The production of photons in high momentum transfer colhsions is useful for a variety 
of reasons [|l|. The single photon invariant cross section is one of several important sources 
of information on the gluon distribution through the role played by the subprocess qg 
0. The photon plus jet cross section may also provide additional constraints as the quantity 
of such data is increased. Photon pair production has already played a role in confirming 
the quark charge assignments and is one of many large momentum transfer processes used 
to check the QCD-based description of high-p^- scattering Most recently, photon pair 
production has received attention because of its role in the search for the Higgs boson in the 
intermediate mass region [Q. The signature provided by the decay i7 77 is particularly 
clean, enabling one, in principle, to detect the relatively rare processes of Higgs production. 
However, a crucial element in this search is to understand the two-photon background that 
results from conventional scattering processes 0. In addition, it has been suggested that the 
observation of the production of photon pairs may serve as a useful signal for new physics 
processes 0. 

The purpose of this paper is to present predictions for the event structure in two photon 
events. In particular, we will focus on events with one additional jet in the final state. 
Such events can be conveniently analyzed in terms of several angular distributions and the 
structure of a two-dimensional Dalitz plot. In addition, the presence of the jet allows one 
to define three classes of events depending on whether the jet has the most, second most, 
or least transverse energy when compared to the two photons. A similar analysis has been 
proposed for events with one photon and two jets in the final state 0]. 

The plan of the paper is as follows. In Sec. |I| a brief description of the program will be 



given. In Sec. |T| the characteristics of the two photon plus jet final states will be presented 
for a set of cuts which could be appropriate for either the CDF or DO experiments. Our 



summary and conclusions will be presented in Sec. IV. 
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II. 0{a'^as) CALCULATION 



The program used for this analysis has been described in Ref . , so only a brief summary 
will be presented here. To order 0{a'^) the production of photon pairs proceeds via qq 77. 
In addition, there are contributions from final states with a single photon and a jet where 
the jet gives rise to a photon. This is estimated by including the qq —>■ 7(7 and qg 75 
subprocesses convoluted with the appropriate photon fragmentation function. The latter is 
formally of 0{a/as) which yields an overall factor of when combined with the factor of 
aag coming from the subprocess. Similarly, one can have final states with two jets, each 
of which fragments to a photon plus hadrons. This contribution is given by the sum of all 
2 — >■ 2 parton-parton scattering subprocesses convoluted with two fragmentation functions. 
These last two contributions are referred to as single and double bremsstrahlung processes, 
respectively. Although formally of one can also have significant contributions from 

the subprocess gg 77 at high energies where small values of x become important since the 
gluon distribution is large enough there to compensate for the suppression from the running 
coupling. 

The 0{a'^as) next-to-leading-logarithm calculation includes contributions from the one- 
loop corrections to qq —>■ '-ff and the 2^3 subprocesses qq — > 77(7 and qg 775. The 
actual calculation is performed using a combination of Monte Carlo and analytic integration 
techniques. The regions of phase space corresponding to soft or collinear singularities are 
isolated by the use of cutoffs, 6s and 6c, applied to the Mandelstam variables for the 2 — 3 
subprocesses. The squared matrix elements are integrated over the singular regions using 
dimensional regularization. The soft singularities cancel with corresponding terms from the 
one-loop corrections. Finally, the collinear singularities are factorized using the MS scheme 
and absorbed into the corresponding distribution or fragmentation function. All remaining 
integrations are performed via Monte Carlo. When properly executed, the cutoff dependence 
in such a calculation cancels between the two-body and three-body contributions which is 
appropriate since the cutoffs merely serve to mark where one switches from a Monte Carlo 
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to an analytic integration technique. Additional details can be found in Ref. [Q. 

The calculation outlined here includes contributions from all 0{a'^as) subprocesses where 
the photons both are part of the hard scatter. There are also fragmentation contributions 
where the hard scattering process is calculated to 0{al) and two of the jets fragment to 
photons. These contributions have not been incorporated into the calculation. However, 
for collider energies isolation cuts are employed to help reduce the hadronic background to 
the photon signals. These isolation cuts greatly reduce the fragmentation contributions. 
Another point to remember is that this calculation is next-to-leading-order for observables 
involving the photon pairs. However, for the two photon plus jet final states only the 2 — >■ 3 
subprocesses contribute, so this represents a leading order calculation for that case. There 
are, for example, corrections from 2 — 4 subprocesses where one of the final jets is not 
detected. These contributions are beyond the scope of this calculation. 

III. PREDICTIONS 

Predictions for two photon and two photon plus jet final states were generated using the 
program described in the previous section. The CTEQ2M distributions P] were used with 
the factorization and renormalization scales chosen to be the maximum Et in each event. 
A minimum transverse momentum cut of 10 GeV is applied to the photons and the jet. In 
addition, the pseudorapidities of the two photons are required to satisfy |?7| < 1 and that of 
the jet to satisfy \ri\ < 3. Both photons are required to be isolated, i.e., to have less than 4 
GeV of additional energy in a cone of radius R = y/Arj'^ + of 0.7 centered on the photon 
direction. Effectively, this just requires the two photons and the jet to be separated by more 
than AR = 0.7. Such a cut greatly reduces contributions from subprocesses involving photon 
fragmentation functions. In the experimental analysis the isolation cut also discriminates 
against jets fragmenting into a leading tc^ or t], and thus greatly decreases the background 
to the true two photon signal. These cuts are typical of what would be appropriate for an 
analysis by either the DO or CDF groups. 
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In Fig. the cross section is shown for two photon and two photon plus jet production 
as a function of the total transverse energy {EjJ) of the two photons. The photon and jet 
(if present) are subject to the cuts described above. The top curve shows the total cross 
section while the bottom curve shows the cross section with a jet satisfying the above cuts. 
The fraction of the cross section coming from events with a jet in the final state that satisfies 
these cuts varies from about 60% for E]? = 100 GeV to about 30% at the upper end of the 
curve. This ratio is larger than might have been expected since the rapidity interval for the 
accompanying jet is much larger than that for the photons. 

In Fig. the fractions of the two photon plus jet cross section coming from the qq and qg 
subprocesses are shown versus the total transverse energy in the final state, denoted by T,Et- 
At low TiEt the qg subprocess dominates while the qq subprocess dominates at high T^Et 
values. This is primarily due to the relative sizes of the gluon and quark distributions at 
low and high values of x. 

For those events with a jet in the final state, the fractions of the cross section are plotted 
in Fig. ^ according to whether the jet has the largest, second largest, or least E^ among 
the two photons and the jet. It is useful to classify the events in this manner since the jet 
is the one object distinguishable form the others. The photon was used in this manner in 
a similar analysis of photon plus two jet events in Ref. 0]. The results in this figure show 
that at high T,Et it is most likely that the jet will have less Ej- than either of the photons, 
while at lower values it is more likely that the jet will have more E^ than one or both of the 
photons. 

For events with two photons and a jet in the final state, four kinematic variables in 
addition to the three-body center-of-mass energy must be specified in order to determine 
the final state configuration. It is convenient to use scaled energy variables for two of these. 
In the parton-parton center-of-mass frame the energies of the final state particles can be 
labelled as E3, E^, and E^ in decreasing order of energy. The energy fractions Xi = 2Ei/Etot 
can then be constructed, with Efot = E3 + E4 + E^. Hence, X3 + X4 + X5 = 2. The variable X3 
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satisfies | < X3 < 1, where the lower hmit corresponds to the symmetric "Mercedes Benz" 
configuration with all three x[s being equal. As X3 tends towards one, the vectors of particle 
4 and 5 become coUinear and opposite to that of particle 3. The allowed range for X4 is 
0.5 < X4 < 1 where the lower limit corresponds to = X5 with particle 4 and 5 being 
coUinear while the upper limit corresponds to X5 = and + X4 = 1. Only two of the x'^s 
are independent; these will be taken to be X3 and X4. Two angles are needed to complete 
the description of the final state. In the parton-parton center-of-mass frame, let 6* be the 
angle between the direction of the incoming proton and the direction of the jet or photon 
with energy E^. The plane containing the beam and this jet or photon will be referred to 
as the scattering plane. The second angle is ip which is the angle between the normal to the 
scattering plane and the normal to the plane containing the remaining jet or photons. 

The main point of this analysis is that examination of various distributions in terms of 
these kinematic variables may provide useful information concerning the accuracy of the 
current QCD-based description of two photon production. For example, some lego plots of 
X3 versus X4 for events where the jet is the most energetic are shown in Fig. ^, |^, and ^. In 
these cases the two photons recoil against the jet. The pr, rj, and isolation cuts described 
above have been applied. In Fig. ^ a cut requiring the total Ej- to be greater than 60 GeV 
has been applied. No obvious structure is visible. However, any structure is hidden because 
the total Ejn cut is low with respect to the individual minimum Ej- requirements. In Fig. 
^ the cut on the total Et has been raised to 200 GeV and in Fig. |] to 400 GeV. As the 
minimum for the total Et is raised an enhancement in the region X3,X4 —>■ 1 begins to 
appear, refiecting the existence of a pole in the matrix element at 1. In this 

region one of the photons becomes soft, i.e., approaches the minimum Et value of 10 GeV. 
For high values of the total Et it is expected that the majority of events in this class would 
have a high Et jet, a slightly less energetic photon, and a relatively soft second photon. 

The lego plots for the case where the jet is the second most energetic show a similar 
structure. A related, but somewhat different, pattern emerges when the jet has less energy 
than either of the photons. For this configuration the jet and the photon with the least 
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energy of the pair recoil against the most energetic photon. Again, a strong pole structure 
becomes increasingly evident for progressively higher values of total Et- However, as shown 
in Fig. 1^, the shape of the emerging pole structure is somewhat broader in this case as 
compared to the case where the jet is the most energetic. This can be ascribed to the fact 
that the subprocess qg 775 does not have a soft singularity when the quark energy goes 
to zero, but it does when one of the photon energies goes to zero. This can be further 
illustrated by comparing Fig. | and Fig. ^ which show the contributions to Fig. |^ from the 
qq — i> ■yyg and qg 775 subprocesses, respectively. 

We have examined the cos(6'*) and ip distributions for the three cases where the jet has 
the most, second most, or least energy with minimum total Ej- cut of 200 GeV applied. The 
ip distributions are essentially the same for all three cases while some minor differences are 
apparent in the cos(^*) distributions. However, the exact shapes are rather cut dependent 
and so have not been shown. The appropriate distributions can easily be generated once 
the appropriate cuts have been specified. 

Several two photon events have been observed by the CDF Collaboration with values of 



total Et exceeding 400 GeV |TD[. Two of the three events have a jet associated with the two 
photons. In one event the jet is the most energetic and for the other it is the least energetic. 
Based on the previous discussion one might expect that the photon or jet with the least 
energy might have an Et value near the minimum allowed, i.e., that the event would be 
in the region of the corresponding lego plot where the enhancement due to the pole occurs. 
Instead, the energies are closer than would have been anticipated. The unusual topologies 
of these events may simply be due to statistical fluctuations, but the continued presence of 
a large number of such events as the data sample grows could signal the need for additional 
sources of high mass photon pairs in the calculation. 
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IV. SUMMARY AND CONCLUSIONS 



As preparations for the next generation of colliders continues, interest in the 77 channel 
has remained high since it can be used to search for the intermediate mass Higgs boson. 
It is, therefore, necessary to understand the conventional QCD production mechanisms for 
photon pairs in order to be able to provide rehable background estimates. The study of the 
event structure of two photon plus jet events using the distributions suggested in this paper 
may well provide useful tests of our understanding. 
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FIGURES 

FIG. 1. The two photon cross section (upper curve) and the two photon plus jet cross section 
(lower curve) versus the total transverse energy of the two photons, Ej^, using the cuts discussed 
in the text. 

FIG. 2. The fractions of the two photon plus jet cross section coming from the 
q'^(solid) and the qg subprocesses (dashed). 

FIG. 3. The fractions of the two photon plus jet cross section where the jet has the highest 
(solid), second highest (dashed), or least (dotted) Et- 

FIG. 4. Lego plot of versus X4 when the jet has the largest of the three Et values and HEt 
is required to be greater than 60 GeV. 

FIG. 5. Lego plot of versus X4 when the jet has the largest of the three Et values and TiEt 
is required to be greater than 200 GeV. 

FIG. 6. Lego plot of X3 versus X4 when the jet has the largest of the three Et values and TiEt 
is required to be greater than 400 GeV. 

FIG. 7. Lego plot of X3 versus X4 when the jet has the least of the three Et values and HEt 
is required to be greater than 400 GeV. 

FIG. 8. The contribution to the plot in the previous figure from the qq initial state. 

FIG. 9. The same as the previous figure except for the qg initial state. 
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